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Abstract 



In this paper we consider the variational approach to cactus trees (Husimi 
trees) and the more common recursive approach, that are in principle equiv- 
alent for finite systems. We discuss in detail the conditions under which 
the two methods are equivalent also in the analysis of infinite (self- similar) 
cactus trees, usually investigated to the purpose of approximating ordinary 
lattice systems. Such issue is hardly ever stated in the literature. We show 
(on significant test models) that the thermodynamic quantities computed by 
the variational method, when they deviates from the exact bulk properties 
of the cactus system, generally provide a better approximation to the be- 
havior of a corresponding ordinary system. Generalizing a property proved 
by Kikuchi, we also show that the numerical algorithm usually employed to 
perform the free energy minimization in the variational approach is always 
convergent. 

KEY WORDS: Cactus tree, Husimi tree, cactus approximation, lattice model, 
Ising model. 



1 Introduction 



Cactus trees are lattices with a branched topology [jl], 0, and usually also 
a self-similar structure Model systems on cactus trees are interesting 

mainly because of two reasons. First they often provide reliable approxima- 
tions to more realistic models on ordinary lattices ^ ^, and second their 
statistical mechanics can be generally worked out exactly ||, If. Be- 

cause of these facts, a lot of physical systems have been investigated in the 
framework of tree lattices: a variety of Ising-like models 0, ||, |10|, |ll 



Potts models |]T2|, spin liquids ||T^, systems with quenched disorder 



polymers |jT5|, |l^, [T^, |T^, |T9|, abelian sandpiles ||2^, electrons in binary al- 
loys 1^ and amorphous solids |2^. The simplest class of lattice models, 
i.e. Ising models, have been most widely investigated also on cactus trees. 
In order to the approximation of systems on ordinary lattices, it has been 
shown by Monroe that the cactus approximation turns out to be particu- 
larly successful in two relevant cases, namely systems with multi-site inter- 
actions 1^ and frustrated systems 0. In both cases the simple mean field 
theory and the Bethe approximation fail in predicting a qualitatively correct 



phase diagram In the special case of the fully frustrated antiferro- 

magnetic Ising model on the triangular lattice, the same holds even for large 
order Cluster Variation Method (CVM) while the cactus approxima- 
tion yields qualitatively correct results Recently Monroe himself has also 
shown that a series of cactus approximations with larger and larger building 
blocks allows not only a more precise determination of phase diagrams ||^ 
but also quite good estimates of critical exponents As a consequence 
of such and other positive results, considerable attention has been devoted 
to the properties of cactus trees and of the methods by which they are 
studied 0. 

In most papers dealing with cactus trees, calculations are based on the 
self-similar structure of the system and this feature is exploited to determine 
its physical properties. In the following we shall refer to such kind of treat- 
ments as to the recursive approach. Nevertheless it is known that cactus 
trees can be studied also by means of a variational approach, equivalent to 
the CVM with a special choice of basic clusters 0. The two methods are in 
principle equivalent and both give the exact solution for finite cactus trees. 
Nevertheless, to the purpose of approximating ordinary lattice systems, one 
is usually interested in determining the bulk properties of an infinite cactus 
tree. This can be done exactly only by means of the recursive approach, by 
evaluating the limit of a recursion relation, or more precisely by investigating 
the attractor of a dynamical system defined by the recursion relation itself. 
In order to employ the variational approach as well, one usually assumes 
some degree of "translational" invariance 0]. Even when such assumption 
is actually verified in the interior (bulk) of the cactus tree, the variational 
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free energy density evaluated in this way turns out to be only the bulk contri- 
bution, not the exact one. The contribution of surface sites, whose number 
increases exponentially as the tree is expanded , cannot be neglected even 
in the thermodynamic limit. This has not always been stated so clearly in 
the literature |0] , and one might expect that a minimization of the bulk free 
energy density yields bulk equilibrium properties. We point out that this is 
true only under certain conditions, namely when the limit of the recursion re- 
lation exists (that is the associated dynamical system has a fixed point, that 
is the bulk invariance condition actually holds [IH), and as far as the equation 
of state is concerned. The latter issue refers to the possibility of multiple 
solutions, i. e. coexistence phenomena, for which the two methods generally 
predict the same solutions (if the fixed point exists) but different first order 
transitions. In other cases (when the fixed point does not exist) the presence 
of boundaries may dramatically affect also the bulk behavior of the cactus 
system, which turns out to be completely different from that predicted by 
the variational approach. By means of actual calculations performed on sig- 
nificant test models we show that, when the variational approach deviates 
from the recursive approach (that is from the exact solution), nevertheless it 
generally yields a better approximation to the behavior of a corresponding 
ordinary lattice system. In the last part of the paper we discuss a particular 
issue related to the variational approach. We show an important property 
of the numerical algorithm generally used to perform the free energy mini- 
mization, known as Natural Iteration Method (NIM). It turns out that the 
free energy decreases at each step. We actually generalize a proof given by 
Kikuchi for the Bethe approximation, coinciding with the variational 
approach to the Cayley tree . 



The paper is organized as follows. In Sec. II we introduce the most impor- 
tant features of cactus trees. Moreover we introduce the variational approach 
and the recursive approach for the case of a finite cactus tree. We show that 
the variational approach is exact because the probability distribution of the 
microscopic state of the system has actually the factorized form predicted 
by the CVM. In Sec. Ill we consider the limit of an infinite cactus tree with 
self-similar structure, taking some restricting hypotheses, that however in- 
clude quite a large number of relevant cases. As previously mentioned, we 
show that, as far as equations of state are concerned, the recursive approach 
reduces to the variational approach, provided the limit of the recursion re- 
lation exists. In Sec. IV we work out the variational approach for three 
test models already investigated by the recursive approach, namely the Ising 
models with pure 4-spin interaction on the square cactus the antiferro- 
magnetic Ising model on the triangle cactus [Q, and the Ising model with 
pure 3-spin interaction on the triangle cactus As previously mentioned, 
we argue that the variational approach provides a better approximation to 
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the phase behavior of corresponding ordinary lattice systems, even if it gen- 
erally does not coincide with the actual bulk behavior of the cactus systems. 
In Sec. V we prove the property of free energy decreasing and discuss some 
consequences of such property concerning the results of the previous section. 
Finally in Sec. VI we give some concluding remarks. 

2 Fundamentals on cactus trees 

We introduce a cactus tree as a lattice with special properties. Let us as- 
sume that each lattice site is characterized by a state variable and that a 
hamiltonian H defines interactions among sites (and of sites with external 
fields) as a function of the state X of the whole system. If interactions have 
finite range, then it is possible to define a set M. of main clusters of sites, so 
that the hamiltonian can be written as a sum over all the main clusters 

H{X) = hMixM), (1) 

where xm denotes the set of site state variables in the main cluster M and 
hui^M) is a part of the hamiltonian depending only on xm- Let us notice 
that the choice of M. (and also of the set of functions {/iAfjMe^vi) Kiay be 
not unique. We require that M C M' can never occur for any M, M' G M.\ 
otherwise we can give another definition of M. to avoid this. Let us also 
define the set of intersections of main clusters 

J = {MnM'{= Jmm') ; M, M' G M}. (2) 

We say that a lattice is a cactus tree if it is possible to define M. non- 
trivially (i.e. not made up of a single element coinciding with the whole 
lattice), so that the following property holds for any integer n > 3 and for 
any Mi,M2, . . . , M„ G 

JmiM2^ • ■ ■ ) JMiMi+1, ■ • ■ 5 JMuMi 7^ =^ 
=^ JmiM2 = . . . = JMiMi+i = . . . = Jm„Mi- (3) 

Eq. (D) states that in the lattice there are no "loops" made up of main 
clusters, and is the reason why the lattice can be referred to as "tree" . It is 
easy to show that such condition implies JCiJ' = for any J, J' G J^. Morita 
refers to clusters in the set JT" as joint dusters , because they actually act 
as "joints" among main clusters, and the intersection of two joint clusters is 
always empty. A possible structure of a cactus tree is depicted in Fig. |l|. 

We now introduce some more definitions and properties, useful to deter- 
mine the factorized form of the system probability distribution (pd) P{X). 
Let us consider two different main clusters M, M'; we say that M and M' are 
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connected if it is possible to find Mi , . . . , M„ G so that Jmmi , • • • , Jm„m' 7^ 
0. We can assume that any M, M' G Ai are connected, otherwise our sys- 
tem would be made up of two or more non-interacting subsystems. Moreover, 
given some joint cluster J, we say that M, M' are connected excluding J if 
there exists a path Mi, . . . , M„ connecting M to M' (in the sense defined 
previously), so that Jmmi, ■ ■ ■ , Jm„m' 7^ J- We shall shortly denote this by 
M ~ M'. It is easy to show that in this way each J defines an equivalence 
relation ~ in which partitions the set of main clusters into equivalence 
classes. It is also possible to prove that such equivalence classes coincide 
with the branches departing from J, defined as 

BjM = {M' eM\M' ^ M}. (4) 

for all main clusters M D J. For example in the system of Fig. |1| J123 and 
<^i45 give rise to three branches, while Jig gives rise to two branches. We now 
define the branch hamiltonians in the following way 

Hjm{xj,Xjm)= Y1 hM'{xM'), (5) 
M'eBjM 

where Xjm denotes the state of the branch minus the base cluster J. We 
also define the partial (branch) partition functions 

Wjm{xj) = E e-'^^-^^^-^"^), (6) 

where the sum runs over all possible states Xjm and (3 = l/ksT, being ks 
the Boltzmann constant and T the absolute temperature. From Eq. 
Eq. @, and the fact that {Bjm}mz)J gives a partition of Ai it turns out 
that, for each joint cluster J, the total hamiltonian can be expressed by 

H{X)= J2 Hjm{xj,Xjm), (7) 

where the sum runs over all main clusters M containing J. Using also Eq. (^, 
this allows to write the pd pj{xj) of the state xj of a joint cluster J as follows 

pj{xj) = n Wjm{xj), (8) 

where Z = Y,x e~^^^^^ is the (total) partition function. Fig. Q shows that 
one more way to write the hamiltonian ([l|) is, for each main cluster M, 

H{X) = hMixM) + Yl HjM'{xj,XjM'), (9) 

JCM M'DJ 
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where the outer sum runs over all joint clusters J contained in M, and the 
inner one over all main clusters M' containing J, except M itself. Using 
again Eq. (|^), the pd Pm{xm) of the state xm of a main cluster M turns out 
to be 

Pm{xm) = Z-'e-"^'"^^''^ II II WjM'ixj). (10) 

JCM M'DJ 

We also notice that 

E i-E(cj-i) = i, (11) 

where cj denotes the connectivity constant of the joint cluster J, i.e. the 
number of main clusters M such that M D J. Such "sum rule" can be easily 
proved by considering any starting main cluster and adding the other ones 
by a "growth" procedure. It turns out that each joint cluster J implies the 
addition of Cj — 1 main cluster, whence Eq. (pAl). 

Taking into account Eqs. (H), ([T0|), (p^, and performing some straight- 
forward manipulations, it is possible to prove that the pd of the whole system 
P{X) = Z^^e'^^^^'^ takes on the factorized form 

n pm{xm) 

Jej 

We can now write the free energy as a function of pds of main and joint 
clusters only. As far as the entropy S is concerned, we split each contribution 
from a joint cluster among main clusters that contain it, and write 

S/kB = {-\0gP{X)) = - J2 (logPM(xM)- E ^^logpj(xj)\, 

MeM \ JcM I 

(13) 

where (■) denotes an ensemble average. For the internal energy f/, from 
Eq. (|I]) we simply have 

U={H{X)) = {hM{xM)). (14) 

Expanding ensemble averages, we can then express the free energy F by 

j3F = j3U - S/kB = E Ypm{xm)^m{xm), (15) 

where 

^m{xm) = PhM^XM) + logpMixM) - bj logpj(xj) (16) 

JCM 
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and 6j = 1 — 1/cj. It turns out that this free energy expression, which is 
exact, coincides with the CVM free energy where the set of basic clusters 



is Ai. Then the exact thermodynamic equihbrium state can be determined 
by the minimization of this free energy with respect to main cluster pds 
Pm{xm), with suitable compatibility constraints for joint cluster pds pj{xj). 
The latter must be obtained as marginal distributions for all M D J by 

Pjixj) = Pm{xm), (17) 

where xm\j denotes the state of the main cluster M minus the joint cluster 
J. This concludes the general discussion as far as the variational approach 
is concerned. 

As far as the recursive approach is concerned, we have to introduce a 
simple relation between the partial partition function of a branch and those 
of its "sub-branches" . Let us consider for instance a joint cluster in Fig. || 
and the hamiltonian of the branch towards the central main cluster. It is 
easy to relate the latter to hamiltonians of branches starting from the other 
two joint clusters displayed. In general it is possible to write 

Hjm{xj,Xjm) = hM{xM) + X! X! Hj>M'{xj',XjfM'), (18) 

J'CM M'dJ' 

where the outer sum runs over all joint clusters J' contained in M except J, 
and the inner one over all main clusters M' containing J', except M itself. 
Hence from Eq. (|^) 

Wjm{xj) = E e-^'^-^^-) n n Wj,M'{x.r). (19) 

^M\J J'CM M'dJ' 

J'^J M't^M 

By means of this equation (in a recursive manner, starting from the bound- 
aries) it is possible to determine partial partition functions Wjm{,xj) for all 
branches Bjm departing from each joint cluster J of a finite cactus tree. 
Eqs. (P) and ([T0|) then provide respectively joint and main cluster pds, from 
which all equilibrium thermodynamic properties can be derived. 

3 The variational and the recursive approach 
to infinite cactus trees 

In this section we consider the limit of infinite cactus tree with self-similar 
structure. This case is relevant to the approximation of model systems on 
ordinary lattices, characterized by translational invariance. We take some 
restricting hypotheses, that however include quite a large number of cases: 
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we assume that (i) site state variables are scalars, and (ii) joint clusters 
are only single sites (we shall simply say joint sites in the following). In 
order to have a self-similar structure, we also require that (iii) each main 
cluster contains the same total number of sites (with the same number n of 
joint sites), and (iv) all main cluster hamiltonians are equivalent. The latter 
condition reads 

hM{x) = h{x) VM G (20) 

with 

X = {xo,Xi, . . . . (21) 

Here x denotes the total state of a main cluster: joint site states are denoted 
by {i = 1, . . . ,n), while Xq denotes the state of the portion of a main 
cluster not covered by joint sites. In principle we assume there are some 
terms in h{x) that make joint sites distinguishable (for instance interactions 
with n different external fields). Defining a set of connectivity constants q 
(i = 1, . . . , n), we build the cactus tree with the usual growth procedure. We 
attach Cj — 1 equivalent main clusters to the i-th joint site of a starting main 
cluster, and then we produce new "generations" iterating the procedure. We 
obtain a structure of the type depicted in Fig. |^, which satisfies all previous 
requirements. Let us notice that such a system turns out to be self-similar 
only in the thermodynamic limit. 

It is not possible to apply to this case the variational approach as de- 
scribed in the previous section, because one would have to deal with an 
infinity of variational parameters. Therefore we assume an invariance condi- 
tion 

p^,{x)=pix) (22) 

for main clusters M in the interior (bulk) of the tree, that is far from the 
surface. In this way the variational approach can only determine the bulk 
equilibrium state, that is p{x), which should approximate that of a corre- 
sponding ordinary lattice model, and is assumed independent of boundary 
conditions. In the hypothesis that Eq. ( P^D holds, in the bulk we shall have 
only a number n of (in principle) different "types" of joint sites, i.e. n differ- 
ent joint site pds Pi{xi) {i = 1, . . . ,n). For convenience we denote the state 
of a bulk main cluster minus the i-th joint site as 

Accordingly joint site pds can be written as marginal distributions by 

As mentioned in the introduction, we only evaluate the bulk free energy 
density / (per main cluster) and perform a minimization of /, expecting to 
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determine the bulk equilibrium state. Taking into account Eqs. (]T5|), (|T6|), 
and the invariance assumption Eq. (p2D , we can write 



^f = J2pi^Mx), (25) 

X 

where 

n 

(^(x) = ph{x) + logp{x) -J2bdogpi{xi), (26) 

i=l 

and bi = 1 — The free energy density / is a functional mp{x) only, being 
Pi{xi) dependent on them via Eq. (^). We then work out the minimization 
with respect to p{x), using the Lagrange multiplier method to impose the 
normalization constraint 

Ep(^) = 1- (27) 



Let us define the functional 



Pfx = /5/ - A 



^p(a;) - 1 



(2^ 



where A is the unknown Lagrange multiplier, to be determined imposing the 
constraint. Taking the derivatives of fx with respect to p{x) (making use of 
Eq. (E5)), and setting them to zero, we obtain 



pix) = z-'e-f"'^''^fl[p,ix,t, (29) 



i=l 



where z is related to A in an irrelevant way. We then simply take the sum of 
both sides of Eq. ( p9|) over all the main cluster states x, and impose Eq. (p7|), 
obtaining 

n 

z = Y^e-^'^^^Y[[p,{x.t- (30) 



Eq. (PPP, with z defined by Eq. (ppj), provides a fixed point equation for 
p{x), which is usually solved via an iterative procedure known as Natural 
Iteration Method (NIM) ||28l. Different solutions may be found starting the 
procedure from different guesses p{x), and the stable phase is determined 
as the solution with the lowest free energy density /. The latter can be 
evaluated by taking the logarithm of both sides of Eq. (pOf), substituting into 
Eq. ( ^61) and then into Eq. (pSf ), yielding 

Pf = -\ogz, (31) 

where z has to be computed at each iteration. In the following we shall 
verify that such a criterion of stability generally does not predict the actual 
first order phase transitions for an infinite cactus system, because it does 
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not take into account surface contributions to the free energy. Nevertheless 
it seems to be the most reasonable way to approximate the phase behavior 
of a corresponding ordinary system 0, ^ . 

Let us turn to the recursive method and consider a branch of our cactus 
tree. It is easy to see that partial partition functions Wjm{xj) depend only 
on the joint site index i and the number k of generations attached to it. For 
a fc-th generation branch, partial partition functions can then be denoted by 



Wi^kixi) (i = 1, . . . ,n), and the recursion relation Eq. (|T^) reads 



Wi,kix,) = ^e-^'^(^) n m',k-iix,)Y''-' . (32) 

i' = l 

Such equation could determine in principle all partial partition functions un- 
til the thermodynamic limit limfc^oo ^i,k{xi)i necessary to obtain bulk prop- 
erties. In practice this is not actually possible, because such limit equals 
infinity. However one usually performs some simple manipulations, leading 
to a feasible recursion relation, but this will be shown in the next section 
for particular examples. By now we only show explicitly that the recursion 
relation Eq. ( p2D is equivalent to the NIM equations (|29| ) and (^) in the ther- 
modynamic limit k —>■ oo, in the hypothesis that such limit exists The 
existence of the limit is actually equivalent to the bulk invariance condition 
Eq. (^21) . In the framework of the recursive approach one usually computes 
the i-th (bulk) joint site pd Pi{xi) by (i) attaching q k-th generation branches 
of the i-th type, (ii) evaluating the central site pd Pi,k{xi), and (iii) taking 
the limit 

Pi{x,i) = lim Pi^kixi). (33) 

fc^oo 

Notice that one actually considers the central site of n (in principle) different 
trees. To be rigorous, this does not evaluate correctly the properties of sites 
close to the surface of the real tree, but becomes exact for bulk sites, still 
in the hypothesis that the limit k ^ oo exists. Specializing Eq. (§) we can 
write 

p,,kix,) = Zr^ [W^,kix,)r , (34) 

where 

Z^,k = E [W^A^^T (35) 

Xi 

is the partition function of the i-th tree made up of q k-th generation 
branches. All ^ tend to the partition function Z of the infinite tree in 



the limit k oo. Substituting Eq. (|34) into Eq. (|3^) , and multiplying both 
sides by [pik{xi)]^' , we obtain 

n n 

P^A^^) = ^i^" n ■Ee'"'^^'^[P.^(^0]'^ mP^',k-lM]'''■ (36) 

i"=l i'=l 
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Let us define 



^ i'=l 



(37) 



Using Eq. (|3|), the fact tliat J2x = Ex^Ex^,, Eq. (|3|), and Eq. (gSj), it is 
possible to prove tliat 

n ^ 
Zi,k = Z^ i^ ' Y\_ 



(38) 



i'=i 



wliicli ensures normalization of Pi,k{xi)- Moreover it is evident from Eq. 
that, if the hmit Eq. (|33|) exists, then 



hm Zi^k = z 

k—^00 



(39) 



(independently of i), where z is defined by Eq. (|30|). As a consequence, 
remembering Eq. (|2^) , we see that Eq. (^6|) in the limit A; cxd is equivalent 
to Eq. (|29| ) "marginalized" to joint site pds (i.e. after a summation of both 
sides over x\i), which proves the equivalence with the NIM. 

By the way we also notice that, taking into account Eqs. (pi]), (PU]), 
and (|38|), we can write 



Pf 



lim 

k—foo 



log Zi^k + 

Ci 



E 

i'=i 



Ci' 



Ci/ 



■ log Zi' 



.k-l 



(40) 



When a Cayley tree is considered, Eq. (^OD reduces to the formula proposed 
by Gujrati (Eq. (3) in Ref. [Q) to evaluate the bulk free energy density in 
the framework of the recursive approach. 



4 Test models 

In this section we consider three test models, which we find very signifi- 
cant and have been previously investigated by the recursive approach P, |[. 
We perform variational calculations and compare the results with those of 
the recursive approach. This should clarify the discussion of the previous 
section and point out analogies and differences between the two methods, 
relating them to both the bulk behavior of infinite cactus trees and that of 
corresponding ordinary lattice systems. 

The first model is an Ising model with pure 4-spin interaction and uniform 
magnetic field on the square cactus (main clusters are square "plaquettes" 
of four sites). Each site is a joint site, characterized by a spin state variable 
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(si, S2, S3, S4 = ±1), and connectivity constants are ci = C2 = C3 = C4 = c. 
The main cluster hamiltonian is 

h{si, S2, S3, S4) = - JS1S2S3S4 - H , (41) 

c 

where J > is the 4-spin couphng constant and H is the magnetic field. 
The NIM equations (^) take the form 

P(S1, S2, S3, S4) = ;2-le-^'^(^^'^^'^3,.4) [P1(S1)P2(S2)P3(S3)P4(S4)]' , (42) 

where p(si, S2, S3, S4) denotes the main cluster pd, 6 = 1 — 1/c, and z is 
determined as usual by normalization. In principle we can distinguish four 
different site pds Pi{si) and four different magnetizations 

rui = {si) = pi{+) - pi{-) (43) 

for i = 1,2, 3, 4, but from the calculation we obtain only homogeneous phases 
with magnetization rrii = m independent of i. 

As far as the recursive method is concerned, on the Ising square cactus 
Eq. (§2]) reads 

W,4si)= E E E e-'''('"''''''''^[W^2,fc-l(s2)W^3,fc-l(53)W^4,fc-l(s4r-' 
S2=±l S3=ltl S4=±l 

(44) 

for site 1, and similarly (by a circular permutation of subscripts) for sites 
2, 3, 4. Eqs. (p^ and (^) as a whole read, for s = ±1 

We can define the ratio 

n = (46) 

independently of i, due to the fact that one assumes homogeneous boundary 
conditions and the main cluster hamiltonian possesses a dihedral symmetry. 
Finally, using Eqs. (^31) , (^) and (|46|) , we can compute the magnetization 
as 

m = lim — , (47) 

where, due to Eq. (|4^), obey the equation 

11 



where a = e'^^^f'^ and d = e^^'^. Eq. ( P8| ) is solved recursively with the free 
boundary condition tq = a, corresponding to a magnetic field equal to H 
acting on all boundary sites. 

In Fig. ^ we report the phase diagrams obtained by both methods for H > 
(field inversion H —H simply implies m — > — m). We have set c = 4, in 
order to approximate the model on the ordinary square lattice. The phase 
diagrams turn out to be qualitatively correct, unlike that obtained by the 
mean field theory (see Ref. ^ for a discussion). We obtain a first order phase 
transition line at if 7^ 0, which separates a phase with lower magnetization 
from a phase with higher magnetization. The line terminates at a critical 
point. According to the previous section, the equation of state is the same for 
both methods, due to the fact that the recursion relation Eq. (|48| ) has always 
a fixed point. On the contrary the phase diagrams are only qualitatively 
equivalent but quantitatively different. This is due to the fact that in the 
variational approach a first order transition is determined by a crossover of 
the bulk free energy densities of two different phases, i.e. two solutions of 
the NIM equations obtained by different guess values. On the contrary, the 
recursive method has a fixed starting point, corresponding to the boundary 
conditions, and detects a first order transition as an abrupt change (driven by 
model parameters) in the attractor of the dynamical system defined by the 
recursion relation. The transition observed in this way is the actual transition 
for the system on the cactus tree. Such difference is hardly ever pointed out 
in the literature §]. As it could be expected, when the two competing 
phases degenerate into one (i.e. at the critical point) the two methods give 
the same result, which confirms the fact that they are equivalent as far as the 
equation of state is concerned. Let us finally notice that the phase transition 
of the cactus system depends on boundary conditions, while that predicted 
by the variational approach is completely insensible to them. This makes the 
latter method not exact in predicting the phase behavior of cactus trees, but 
more suitable to approximate (translationally invariant) systems on ordinary 
lattices, which is very well verified in the present case. Let us notice that 
the Ising model with pure 4-spin interaction on the ordinary square lattice is 
self-dual [^, and phase transitions should occur (if they do) along the line 
given by 

sinh(2/3J) sinh{2f3H) = 1. (49) 

From Fig. |^ we observe that the transition line obtained by the variational 
approach nearly coincides (within a tolerance ~ 10~^) with the self-dual line. 
In order to obtain better approximations to ordinary lattice systems also in 
the framework of the recursive approach, an alternative criterion has been 



proposed by Monroe for the location of first order transitions. Such 
criterion is based on the evaluation of the derivative of the recursion relation 
at its fixed points, and has been verified to give the same numerical result 
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as the variational approach, in the case of the Potts model. Based on the 
results of Ref. |^ we can state that the equivalence is numerically verified 
also for the present model, because the transition line obtained by Monroe's 
criterion nearly coincides with the self-dual line as well. Nevertheless an 
analytic proof of the equivalence of such criterion with the minimization of 
the bulk free energy density has not been given yet. As previously mentioned, 
another criterion, based on the evaluation of the bulk free energy density has 
been proposed by Gujrati in the form of Eq. (^Ol). Unfortunately such 



an expression is quite difficult to evaluate numerically in the framework of 
the recursive approach, because it involves a difference between quantities 



tending to infinity in the thermodynamic limit. On the contrary Eq. (|3lD 
shows that the bulk free energy density comes out in a natural way from the 
variational approach. 

The second test model we consider is the antiferromagnetic Ising model 
with uniform magnetic field on the triangle cactus (main clusters given by 
three site plaquettes). Each site is a joint site, characterized by a spin state 
variable (si, S2, S3 = ±1), and connectivity constants are Ci = C2 = C3 = c = 
3. The main cluster hamiltonian reads 

h{si, S2, S3) = - J (S1S2 + S2S3 + S3S1) - 7/ fllL^llL^^ (50) 

where J < is the antiferromagnetic coupling constant and H is the mag- 
netic field. The NIM equations are similar to Eq. (^2]) 



p{su S2, ss) = z-ie-^'^^^^'^^'^^) [pi{sMs2)p3{s3)t , (51) 

(with obvious meaning of symbols), while magnetizations can be obtained 
by Eq. (^3]). From the calculation we obtain a homogeneous phase and a 
symmetry-broken phase, where on every triangular plaquette we have (for 
H > 0) two sites with equal positive magnetization and one site with nega- 
tive magnetization. The situation is inverted for H < 0. The phase diagram 
is displayed in Fig. ^ and is symmetric with respect to H = 0. The transi- 
tion line is always first order. This model turns out to be interesting as an 
approximation of the antiferromagnetic Ising model on the ordinary trian- 
gular lattice, for which, due to frustration, ordinary mean field like approx- 
imations [^, included the CVM |2S|, fail in predicting the (qualitatively) 
correct phase diagram, and show a phase transition at zero field. 
As far as the recursive method is concerned, Eq. (|32|) reads 

WlAsi)= E E e"^'^'"'''''MW^2,fc~l(s2)W^3,.-l(^53)]'='' (52) 

S2=±l Ss=±l 

for site 1, and similarly (by a circular permutation of subscripts) for sites 
2, 3. The procedure is analogous to the previous case, except the fact that we 
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preserve the dependence on i, in order to be able to consider inhomogeneous 
boundary conditions (if homogeneous boundary conditions are imposed the 
dependence on i disappears because of the dihedral symmetry of the main 
cluster hamiltonian) . Being 

"''"'w^ki-y ^^^^ 

we obtain the following recursion relation 

rik = CL 7^ , (54) 

for site 1, and similar ones for sites 2,3 (circular permutation). Magnetiza- 
tions are computed by 

rrii = lim — -. (55) 

The results of the recursive approach turns out to be dramatically affected 
by boundary conditions for the present model. We consider a fixed temper- 
ature UbT /\J\ = 1 and vary the field H, considering the following cases. For 
uniform free boundary condition ri^o = ''^2,0 = '"s.o = « (magnetic field equal 
to H on all boundary sites) we obtain the results displayed in Fig. ^(a). 
The dependence on i is removed but, in a region H G (0, He), the recursion 
relation has no fixed point and displays a limit cycle of period 2. The mag- 
netization of the central site oscillates between the two values shown in the 
figure, the positive value for even k and the negative one for odd k. In both 
cases triangular plaquettes of consecutive generations alternatively display 
two sites with positive magnetization and one site with negative magnetiza- 
tion, or vice-versa. On the contrary for if = and H > He a. fixed point 
exists and a paramagnetic phase with uniform magnetization is obtained. 
The latter is equivalent to that predicted by the variational approach. We 
also consider the case of inhomogeneous boundary conditions ri^ = r2,o = ^ 
and rs^o = (magnetic field equal to if on 2/3 of boundary sites and —H 
on the remaining ones). We obtain the results displayed in Fig. ^(b). The 
behavior is equivalent to the previous one, except in a subinterval of (0, He), 
where a fixed point exists, and the dependence on i is preserved. More 
precisely we obtain the same symmetry-broken phase predicted by the varia- 
tional method, with the same numerical values of magnetizations. This is in 
agreement with the discussion performed in the previous section. We finally 
compare the above results with those obtained by Monroe |^ by solving the 
recursion relation Eq. (|5^ ) in a "sequential" way. Even if this actually corre- 
spond to a sligthly different system (a tree with a "ragged" surface), it turns 
out that a fixed point always exists, and the behavior of magnetizations. 
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displayed in Fig. |^(c), is quantitatively equivalent to that predicted by the 
variational approach (except for the positions of transitions between the uni- 
form phase and the symmetry-broken phase). After the discussion of these 
results, we conclude that there are some cases in which the bulk behavior of 
cactus systems does not provide reliable information about the behavior of 
the corresponding ordinary lattice system, hence it must be employed with 
some caution. Apparently contradicting ourselves, we have to remark that 
the "sequential" recursive procedure worked out by Monroe gives, with re- 
spect to the "normal" recursive method, a transition line which is closer to 



the Monte Carlo result [^] for the model on the ordinary triangular lattice 
(see Fig. Nevertheless this seems to be a peculiarity of the present model. 

The third test model we investigate suggests that the recursive method 
alone may lead to incorrect conclusions about the physics of the model on 
ordinary lattice, if no other information is available. We consider an Ising- 
like model with pure three-spin interaction and uniform magnetic field on 
the triangle cactus (introduced above). This model has been previously in- 
vestigated by the recursive method 0, with the aim of approximating a 
model with three-spin interaction on upward pointing (or downward point- 
ing) triangles of an ordinary triangular lattice. The main cluster hamiltonian 
reads 

h{si,S2,S3) = -JS1S2S3 - H , (56) 

c 

where J > is the plaquette interaction, H is the magnetic field, and c = 3 . 
All calculations are analogous to the previous model. The recursion relation 
turns out to be 

^^^^ = ^ 9 c-i c-i TT^i 7^ 7 ^^^^ 

for site 1, and similar ones can be derived (by the usual circular permutation) 
for sites 2, 3. We obtain the phase diagram shown in Fig. ^ For H > only 
homogeneous phases are obtained, with a first order transition line which 
separates a phase with lower magnetization from a phase with higher mag- 
netization. The line terminates at a critical point. The phase behavior is 
qualitatively similar to the square cactus model described previously. Sim- 
ilarly the two methods predict different first order transition lines but the 
same critical point. On the contrary for H < a. symmetry-broken phase 
appears. According to the variational method, each triangular plaquette has 
two sites with equal negative magnetization and a site with positive magneti- 
zation. This phase is separated from the paramagnetic phase by a first order 
transition line. In almost the same region of the phase diagram the recursive 
method displays a peculiar behavior, involving limit cycles of high order and 
transitions to chaos [Q. In Fig. |^ we report the boundaries of such region. 
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drawn from data published in Ref. 0. We conjecture that in this case the 
correct phase diagram of the ordinary lattice model (3-spin interaction on 
upward/downward triangles) is that predicted by the variational approach 
(cactus approximation), while the anomalous behavior observed in Ref. |^ 
(and also Ref. ||T0| for a similar model) is a peculiarity of the cactus tree. Our 
conjecture is also supported by the fact that, applying Monroe's sequential 
procedure to this case (having some analogies with the previous one), we 
have obtained results in agreement with the variational approach. 



5 Convergence of the NIM 

In this section we examine an important property of the NIM equations (pQ]). 
As previously mentioned, the NIM is a numerical iterative minimization of 
the variational free energy density. By a generalization of the proof given 
by Kikuchi for the Bethe lattice p8[, it turns out that the free energy de- 
creases at each iteration, and the algorithm is always convergent. Let us give 
the proof and then discuss some consequences. Starting from Eq. (pSj), we 
write the difference between the free energies of two consecutive steps of the 
iterative procedure as 

Pif - /) = E iPi^Mx) - p{xMx)] , (58) 

X 

where a hat denotes the latter step, and accordingly 

n 

0{x) = (3h{x) + logp(x) - ^hi logpi(xi), (59) 

i=l 

while '^{x) is defined by Eq. (^Bj). Taking the the logarithm of both sides of 
Eq. (|29|) (where the left side is now denoted by a hat), we can write the NIM 
equations in two different ways, that are 

n 

\ogp{x) = - \ogz - Ph{x) + ^bi\ogpi{xi) (60) 

1=1 

and 

n 

^bi\ogpi{xi) = \ogz + (3h{x) +\ogp{x). (61) 

i=l 

We substitute the former into 0{x), the latter into (f{x), and finally both 
into Eq. (^), yielding 

Pif - /) = E pi^) log + E E M^^) log fl^- (62) 
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Let us now consider the inequahty log^ < ^ — 1, that holds for all real 
numbers ^, and observe that the equality holds only if ^ = 1. By applying 
this argument to both logarithms in Eq. (^) and remembering that all pds 
are normalized at each step, we can finally write 

/-/ < (63) 
/ — / = <^=^ p{x) = p{x) Wx. (64) 

Eq. (p3D means that the free energy can be decreasing or constant during the 
procedure, while Eq. (0) assures that it is constant only if the procedure has 
already reached convergence (i.e. the free energy can only decrease during 
the procedure). 

The above property, generally desirable for a numerical method that aims 
to minimize a function, has some relevant consequences and makes a signif- 
icant difference between the NIM and the recursive method. It is evident 



that Eqs. ( |63D and (p^ ) prevent the dynamical system defined by the NIM 
equations from having limit cycles. In this way the variational approach al- 
ways determines the best solution compatible with the invariance condition 
Eq. (|2^) , but cannot detect whether such hypothesis is too restrictive or 
not. The variational approach describes correctly a symmetry breaking with 
a period less than or equal to the width of a main cluster, but would not 
be able to indicate the existence of phases with higher periodicity (or even 
incommensurate phases, as observed for instance in the ANNNI model [0). 
We then conclude that also the variational approach must be used with some 
caution in the approximation of ordinary lattice systems, if there are reasons 
to suspect that a violation of Eq. ( ^21) occurs, not only in the cactus system. 



In such cases the recursive approach is essential []TT| because, as we pointed 
out in the previous section, the dynamical system defined by the recursion 
relation is able, through the features of its attractor, to indicate the nature 
of the symmetry breaking. 



6 Conclusions 

In this paper we have discussed several properties of the variational approach 
to cactus trees, performing comparisons with the more usually employed re- 
cursive approach. First of all we have pointed out that the variational ap- 
proach is based on an exact factorization of the state probability distribution, 
and can in principle solve exactly finite cactus trees, as well as the recur- 
sive method. Moreover we have considered different issues, concerning the 
bulk behavior of infinite (self-similar) cactus trees and the approximation of 
ordinary lattice systems. 

We have shown that the variational method allows a simple evaluation 
of the bulk free energy density. The minimization of bulk free energy yields 
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the correct equation of state for the interior of the cactus tree. In presence 
of multiple solutions, i.e. coexistence phenomena, the first order transitions 
determined by the variational method are not the exact ones for the cactus 
trees, but turn out to be independent of boundary conditions, and pro- 
vide a reliable approximation (cactus approximation) to phase transitions 
of corresponding ordinary systems. On the contrary the recursive method 
determines the exact bulk behavior of infinite cactus trees, on the basis of a 
change in the attractor of a dynamical system, driven by model parameters. 
Unfortunately such behavior is strongly dependent on boundary conditions, 
and usually provides poorer approximations to ordinary systems, or even 
incorrect results. 

We have also considered the convergence property of the algorithm gener- 
ally used to perform the free energy minimization in the variational approach. 
From the point of view of numerical analysis, we deal with a simple iterative 
solution of a set of fixed point equations, actually the same kind of problem 
which is solved in the framework of the recursive approach. Nevertheless 
the peculiar form of the variational equations allows to prove that the free 
energy corresponding to the current thermodynamic state decreases at each 
iteration and remains constant only if the algorithm has already reached 
convergence. We have shown on a counterexample that this is not always 
true for the recursive method. Of course the free energy decrease is a nice 
property in most cases, because we are interested in free energy minima. 
Nevertheless, on the basis of results obtained on a test model, we have sug- 
gested that it may also lead to missing some complex physics, coming out in 
the form of limit cycles and/or chaotic attractors in the recursive approach. 

In the whole paper, sometimes only by recollecting already known re- 
sults, we have tried to give a unified picture of the two different approaches, 
pointing out analogies and differences. 
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Figure 1: Possible structure of a cactus tree: Mj denote main clusters, Jjj... 
denote joint clusters. The notation is such that J;y...fc = Mi fl Mj — . . . — 
Mk n Mj. The central cluster Mi contains J123, J145, and Jiq. 



Figure 2: An example of (planar) cactus tree with equivalent main clusters, 
and n = 3 different types of joint sites. Connectivity constants are Ci = 3 
and C2 = C3 = 2. Shaded areas denote portions of main clusters not covered 
by joint sites. 
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Figure 3: Phase diagram of the 4-spin Ising model on the square cactus 
(temperature vs. magnetic field). A dashed line denotes the first order tran- 
sition, computed by the recursive method with free boundary conditions. A 
circle denotes the critical point. A thick solid line denotes the same tran- 
sition evaluated by the variational method. A thin solid line represents the 
self-dual line. 



Figure 4: Phase transitions of the antiferromagnetic Ising model on trian- 
gle cacti (temperature vs. magnetic field). The symmetry-broken phase is 
denoted by -|- H — . A thick solid line denotes the first order transition to 
the uniform phase, as predicted by the variational method. A dashed line 
denotes the same transition obtained by the "sequential" recursive method 
(see the text). Squares denote results from Monte Carlo simulations of the 
model on the ordinary triangular lattice (the thin solid line is an eyeguide). 
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Figure 5: Magnetizations of the antiferromagnetic Ising model on triangle 
cacti at fixed temperature A;bT/| J| = 1 as a function of the magnetic field. 
Results obtained by the recursive method with: (a) free boundary conditions; 
(b) reversed field on 1/3 of boundary sites; (c) free "ragged" boundary condi- 
tions ("sequential" method). Solid lines refer to fixed point magnetizations, 
dashed lines to limit cycles (period 2). 
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Figure 6: Phase transitions of the 3-spin Ising model on the triangle cactus 
(temperature vs. magnetic field). Thick solid lines denote first order tran- 
sitions computed by the variational method. A circle denotes the critical 
point. A thin solid line represents the self-dual line. A thick dashed line 
denotes the first order transition between homogeneous phases obtained by 
the recursive method with free boundary conditions. The symmetry-broken 
phase region is denoted by H . Squares mark the boundary of the re- 
gion in which the recursive method displays limit cycles and chaos (the thin 
dashed line is an eyeguide). 
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